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Abstract 

The analytical values of the sunrise master amplitudes were found for one mass 
equal to zero, two other arbitrary masses and an arbitrary external four momentum. 
Differential equations in the square of the external four momentum and masses were 
used to obtain the master integrals. 
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1 Introduction. 



High precision measurements in elementary particle physics require more and more precise 
calculations of multi-loop Feynman diagrams allowing in the same time detailed tests of 
the existing models. 

The widely used integration by part technique allows for reduction of the compli- 
cated task of calculation of radiative corrections to calculation of a limited number of 
scalar integrals called master integrals (MI). Whenever possible it is of great importance 
to know analytically that basic building blocks, as it safes enormously the CPU time 
required for their numerical evaluation. Investigations of the sunrise two-loop diagram 
started a long time ago and it is beyond the scope of that paper to present the complete 
literature of the subject, so we will give here only a partial guideline of the literature 
sending the reader to the references of the quoted paper for further studies. 

The MI of the sunrise graph were found to be expressible as a combination of Lauricella 
functions in H . Their analytical expansions at p 2 = and p 2 = oo were given in || and 
[||. Their values at threshold p 2 = — (mi + m 2 + m 3 ) 2 and pseudothresholds (p 2 = 
— (mi — m 2 + m 3 ) 2 , p 2 = —(mi + m 2 — m 3 ) 2 , p 2 = —(—mi +m 2 + m 3 ) 2 ) were presented in 
U|]. The analytical expansions at pseudothresholds were found in Q. A semi-analytical 
expansion at threshold was given in 0] and also using the configuration space technique 
in |7|, while the complete analytical expansion at threshold was presented in H. Some 



special values of the sunrise MI were obtained in M and [10 . 



2 Calculation of the sunrise master integrals. 

The scalar integrals related to the sunrise two-loop graph can be defined as ( 



A(n,m 2 1 ,ml,m 2 3 ,p 2 ,-ai,-a2 J -a 3 ,Pi,p 2 ) = ^ 2 tt)"- 2 ) 2 



d n ki J d n k 2 



(p-kiY x (p-k 2 f 2 



(k\ + m 2 ) Q i (k\ + m 2 2 Y 2 ((p - h - k 2 ) 2 + mj)^ ' 

(1) 



where (i = 1,2,3 ) are the masses associated with internal lines, p is the external 
momentum, ki,k 2 are loop momenta and ctj, (i = 1,2,3), j3j, (j = 1,2) are integer 
numbers. The integrals are to be performed in n- dimensional Euclidean space. It means 
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we use the dimensional regularization and have performed the Wick rotation. The scale 
parameter \i associated with the dimensional regularization was set to 1. Final results 
can be easily rewritten in Minkowski space by changing p 2 — ► — p 2 . 

The four independent Master Integrals (MI) were chosen as 



rp I 2 2 2 2\ 

F j (n,m 1 ,m 2 ,m 3 ,p 



((2< 



n-2^2 



fl n h ( d n h 0\ 

1 J 2 (kj + m 2 Y^){k 2 + m 2 Y^){{p -h- k 2 ) 2 + m|)«3(i) ' {Z) 

with j = 0, 1, 2, 3; i = 1, 2, 3; a;(0) = 1; a { (j) = 1, for j ^ i; a^j) = 2, for j = i. While 
their (n — 4) expansions are defined by 



I7i / 2 2 2 2\ 1 rp(— 2)/ 2 2 2 2\ , ^ rri(— 1) / 2 2 2 2\ 

Fjin.m^m^m^p ) = _ 4 )2 i j ( m n m 2» m 3»P ) + ^ _ 4 ) F i K^i^iPj 

+ F, (0) (m?,m2,m2,p 2 ) + ... , j = 0,l,2,3 (3) 

The calculation of the MI for m3 = presented here make use of two systems of 
differential equations satisfied by the MI of the sunrise diagram given in Q. One of them 
in the variable p 2 



p2 ^p F ° = ( n - 3 ) F o + m i F i + m l F ^ + m l F 3 

d 3 

p 2 D(m 2 ,m 2 ,m 2 ,p 2 ) — F t , = ^M^ + T, , «=1,2,3 , (4) 

dp 2 po 

where the explicit form of functions T\ (expressed as functions of T ) and My (polynomials 
of p 2 , m 2 , m|, m| ) can be found in ||. The function D is defined by 



D(m 2 ,m 2 ,m 2 ,p 2 ) = (p 2 + (mi + m 2 + m 3 ) 2 ) [p 2 + (m x + m 2 - m 3 ) 2 ) 

• (p 2 + (m x - m 2 + m 3 ) 2 ) (p 2 + (m x — m 2 — m 3 ) 2 ) , (5) 



and the aforementioned function T by 
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T = / h? + = (n-2)(n-4) C(n) ' (6) 



where 



C(n) = (2v^) (4 n) T (^3 - ^) and C(4) = 1 . (7) 

The second system of differential equations, with /Tig cLS cL variable, can be easily 
obtained using expressions presented in the Appendix of [[|. They represent 

A(n, , ml, m ii J° 2 i — 3, —1, —1, 0, 0) and A(n, ml, m\, m\,p 2 , —1, —2, —2, 0, 0) as linear 

combinations of the MI. Remembering that 



A(n, ml ml, m 2 3 ,p 2 , -1, -2, -2, 0, 0) 
A(n, m\,m\, m\, p 2 , -1,-1, -3,0,0) 



<9F 2 (n, mf , m 2 ,, m\,p 2 ) 
dm 2 

1 dF^{n, m\, m\, m 2 ,p 2 ) 

2 <9mo 



etc. 



they are also differential equations in the variable m\ (or in any other mass by means of 
permutation of masses). They can be written as 



8 

tti / 2 2 2 2\ a ri / 2 2 2 2\ , /i tti / 2 2 2 2\ 

F (n,m 1 ,m 2 ,m 3 ,p ) = A 0fi F (n,m 1 ,m 2 ,m 3 ,p ) + A 0jl F 1 (n,m 1 ,m 2 ,m 3 ,p ) 



(^777,3 



+A , 2 F 2 (r7, m 2 , m 2 , m], J7 2 ) + ^3^(77, m?, m 2 , m 2 ,p 2 ) + B 





;Fi(n,m 2 ,ml,m 2 3 ,p 2 ) = A ifi F (n,ml,m 2 2 ,m 2 3 ,p 2 ) + A iA F x {n,m\,m\,m\,p 2 ) 



dm 2 

+A i)2 F 2 {n,ml,m 2 2 ,m 2 3 ,p 2 ) + A ii3 F 3 (n, ml, m\, m\, p 2 ) + B, 



(9) 



where % = 1,2,3 . We do not present here for short the explicit form of the coefficients 
Ai t j and Bi as they can be easily found from [[|. 
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To obtain the values of the MI at m\ = we write the expansion of Fa(n, mf, m\, m\,p 2 ) 
at that point (the expansions of the other master integrals are obtained using Eq.(|8])) 



F^mlmlmlp 2 ) = £ (m^ -J]W > ( 10 ) 

«,;GA fc=0 

where A is a set of all allowed values of a*. The allowed values of ctj one can easily get 
putting Eq . (|ToD and its mf derivatives into Eq.(|]). We have found that there exist only 
two allowed values of the powers a;j and the expansion of F (n, mf, m|, m|,p 2 ) reads 



T7T / 2 2 2 2\ 

F (n,m 1 ,m 2 ,m 3 ,p ) 



oo 



fc=0 



fc=0 



11' 



It consists of two independent series, the regular one denoted as r and the singular 
one denoted as s. It means that the m\ 
equations. The equations Eq.( 



point is a regular singular point [11] of that 
provide also with the following differential equations 



C 2 (n)(m 2 )( n ~ 4 )/ 2 (m 2 + m 2 + p 2 ) 



+ 



2(n-2)(ra-4) 



2 



R 2 ( m lml-p 2 )JL F s = {n _ 3)(m 2 _ m 2 + p2)F s _ 

C 2 (n)(m 2 )( n - 4 )/ 2 (m 2 + m 2 + p 2 ) 
+ 2(n-2)(n-4) 2 ' (12) 



where 



R 2 (—p 2 , mf, m 2 ) = p 4 + m 4 + + 2m 2 p 2 + 1m\p 2 — 2mlm\ . (13) 

On the first sight it is not necessary to know the singular part when looking for the 
value of the MI at m% = as it vanishes in that limit. However as it will be clear from 
the discussion below knowledge of its poles in (n — 4) provides an essential simplification 
in the calculation of the regular part. 



4 



The equations Eg . (|T2|) can in principle be solved, but it is much easier to use system 
of eguations Eq.(f|) with p 2 as an independent variable for calculation of Fq. Substituting 
F (n, m\, m\, m\,p 2 ) and other MI with their mass expansions in the system of eguations 
Eg.(^) we have found that Fq satisfies the following differential eguation 



d 1 1 

dp 2 ' 2p 2V } p 2 R 2 {mi,ml-p 2 ) 



[n - 3) (p 2 {m\ + ml) + (m 2 - m 2 ) 2 ) F { 



+ 2 (, i ^-4) 2 '"^ , "' 2)/2 '" 2 -'"- + ^> 



(14) 



We have found a solution of that differential eguation in a form expanded around 
n = 4. Only pole parts of the function were calculated as only the poles are necessary for 
our purpose of calculating the values of the MI at m^ = 0. The solution reads 



C 2 {n) 



+ 



C 2 {n) 



4(n-4) 2 16(n-4) 

' r>/ 2 2 2 

R(m 1 ,m 2 , -p _ 



2 2 / 2 \ 

6 - \og{m{) - log(m^) H L — — - log ' 



(n-4) 



C 2 (n) St 
-—- log (t) - — 
8p 2 p l . 



p- 



+ 0((n-4) ( 



mt, 



(15) 



where S\ is an integration constant. The unknown function of m\ and m 2 , which remains 
after p 2 integration, was reduced to a single numerical constant Si using the eguations 
Eq-flTJ). 

The system of eguations Eg.(^|) gives no information about Fq, provides however an 
expression for the second coefficient in the regular series. It can be expressed as a linear 
combination of Fq and its derivatives 



F{ 



1 



R 2 (m 2 , m\, —p 2 
4 



[n 



4) 



(n - -)(n- 3) (ma + m\ +p 2 )Fq 



(n-4) 
1 



m\(n — 3)(p + m 1 ) 



JL F r 

dm 2 



[n 



4) 



m\{n — 3)(p 2 + m\) 



JL F r 

dm 20 



+ 



(n - 4) : 



■(m?)("- 2 )/ 2 (m2)(- 2 )/ 2 C 2 (n) 



(16) 
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Putting the expansion Eq.fllTD into the equations Eq.(|4]) one finds a system of three 



differential equations satisfied by Fq, q^iFq and 



dp 2 

±J_ F r 

dp 2 dm 2 



d 8 
dp 2 dm; 



■F n r 



+ 



+ 



+ 



+ 



+ 



+ 



pZ pZ pZ jj m z 



d m\ 








FX 



1 



R 2 (m 2 , m 2 ,, — p 2 
dm 2 \ p 2 



p 2 + m 2 



2p 2 

1m\m\ 3mf 
p 2 2p 2 
2m\ 

p2 

2 2 

m\m2 



p 2 dm 2 
4 + 7(n-4) + 3(n-4) 2 )F c 



'o 



4 2 

nij p 



P 

ml 



2m? 



7m 



+ - 3m? 



6 2 v p* 

C 2 (n)( m 2)(«-4)/2( m 2)(n-2)/2(_ p 2 + m 2 - m 



4(n - 4) 2 p 2 



1 



R 2 (m 2 , m 2 ,, — p 2 
9 



p- 



mj + m 2 



2p 2 



(n-4) 




-F r 

Dm 2 



2m 2 m; 



5m 2 



+ (n-3) 



2 v P 
3m 2 m 2 

2p 2 
9m 2 



p z 

2 2 

m^m^ 



m\ 
' V 
2m^ 

P 2 



+ -i - 3m 2 



4 + 7(n-4) + 3(n-4) 2 F ( 



3m| p 2 2 

!T H "^9 

2p 2 2 2 



m. 



2m 2 



P 

,2)(n-4)/2, 

4(n — 4) 2 p 2 



C 2 (n)(m 2 )(™- 2 )/ 2 (m 2 )( n - 4 )/ 2 (-p 2 - m 2 + m 2 ) 



(17) 



Trying to find Fq from Eq.(|T7j) only would mean to solve a third order differential 
equation satisfied by Fq. That task can be simplified enormously by careful examination 
of the Eq . (|16|) . One finds from it that in the (n — 4) expansion of F[ the terms ~ (n — 4)° 
of Fq contribute to pole terms of F{. As a result, knowing the pole terms of the singular 
part of the expansion (Eq.(|i~5|)) and the exact pole terms of the complete master integral 
F ||, we can find an additional relation between Fq and its mass derivatives. Defining 
the (n — 4) expansion of the Fq by 
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in 



4) 



2 r 0, 



+ 



[n 



4 )V 



+ itf + 0{n - 4) 



(18) 



one finds 



= (ml + m 2 2 + p 2 )F^ - m\(j? + m\)-^F r ^ - m 2 2 (p 2 + m\)^F^ 



+R(ml ml -p 2 ) ( log(t) - ^ (p 2 + 2m? + 2m 2 + - 2 (m\ - m 2 " 2 

+C 2 (n) ( ^ (l3 - 2 log(m?) - 2 log(m 2 )) + ^ (41(m? + m 2 ) 

— log(m 2 )(14m 2 + 6m|) — log(m 2 )(6m 2 + lAm\) \ H -(m\(?>m\ — ml) 

J 6 1/r V 

)/ \ 7Tl 2 7Tl 2 { ^ 2 

(log(m 2 ) - log(m 2 )J + — ^Uog(mi) + log(m 2 ) 

+ — (12m 2 m 2 + 7m? + 7rr4) - log ^^ (12m 2 m 2 + 5m? + 3mj) 
32 64 

l0g(m ' } (12m 2 m 2 + 3m? + 5m^) ) (19) 



64 
where 



p 2 + (mi + m 2 ) 2 - ^/p 2 + (mi - m 2 ) 2 ^ 



p 2 + (mi + m 2 ) 2 + \jp 2 + (mi — m 2 ) 2 

It allows us to reduce the third order equation obtained from Eq.flTTD to a relatively 
simple second order differential equation 



d 2 



d{p 2 ) 



2\2 A 0,0 



p 2 dp 2 °'° 
C 2 (n) 



(p 2 ) 



TDl 2 2 2 

: R(m 1 , m 2 , -p 



'C 2 (n) 



16 



log(t) 



-5, 



32(p 2 )^ 



m 2 — m 2 ,) log — | I — k^m 2 ,)^ 2 — log(m 2 )p 2 H 



m$ 



.(21) 



Its integration is elementary giving dilogarithms as the most 'complicated' functions 
and the solution after all integration constants were fixed reads 
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+ 
+ 
+ 



F o (o) K,m^,0,^)=F ( 



0,0 



C\n) 



m 2 m 2 log {t) 
16p 2 



+ 



— (-m\m\ log 2 (t\) + K - mj) log (t 2 



P 2 A / 2x,, / 2n 13\ log(t) /m? \ 
+ ^^ 2 )+log(mO- T J-3^ r ^y(- + m 1 m 2 j 



log(t) /mj 
32(t-t 2 ) Vto! 
log 2 (t) 



+ mim 2 



32 



.R(m 2 , m 2 , —jo 2 ) + 2m 2 + 2m 2 



32 



-(m 2 + mj) + — 



2 T^(Li 2 (l-f)-Li 2 (l-f; 
16 V ti £ 2 



2 log(ti) (m 2 — m|) + m 2 + m 2 . 



32 



log 



(ti-t)(t 2 -t) 

p2 



log(m 2 ) 
log 2 {ml) 



128 



(13m 2 + mj) + 
(5m| + m 2 ) — 



log(ml) 2 2 A log 2 (m 2 ) 



-(13m 2 + mj) - 



64 

log(m 2 ) log(m2 
64 



128 



(5m 2 + m 2 ,) 



+ 32 K m i + m 2J 



(22) 



where 



mi m 2 
ti = — and t 2 = — . 

m 2 mi 

The constant Si was fixed using the analytically known || result for F (n, 0, 0, m|, p 2 ), 
while two unknown functions of mi and m 2 coming from the integration of the second 
order equation Eq.(^TJ) where fixed from known || results for _Fo(n, m 2 , m|, m\, 0) and 
^F (n, m\, m|, m\,p 2 ){p 2 = 0). The value of p 2 is to be understood, whenever necessary, 
as p 2 — ze, with e an infinitesimally small positive constant. 

Having F^\m\, ml, 0,p 2 ) one can easily get other two independent MI: Ff\m\, ml, 0,p 2 ) 
and F2°\m1, ml, 0,p 2 ). We present here only one of them as the other one can be found 
by an exchange of the masses mi and m 2 . It reads 



Fi (0) (m 2 ,m 2 ,0,p 2 ) 



C 2 



n 



1 

32 + 



R(m\, ml, -p 2 ) log(t) log(t 2 ) 



16p 2 



32p 2 



ml - f (log(t 2 ) + 1) 
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+ 



log(t) 

32 
log(m; 



2 + 



logg) 

p2 



128 



(p 2 + 2m 2 )) - ^| (log(m 2 ) + log(m 2 )) 

5 log(m 2 ) 



+ 



1 

32 



2 log(m 2 ) + log(m 2 ) J + log(m 2 ) ( ~ 
(h-t)(t 2 -t) 



128 



(2-^ + log(tf))log 



Li 2 (l - 



(23) 



We have tested numerically the results of Eq. (|22|) and Eq. (p3|) against results obtained 
by newly developed program jl2], which solves numerically the system of differential 
equations Eq.(f|). In the program one cannot put however one of the masses to zero, 
though very small values of the masses are allowed. Putting sufficiently small masses to 
eliminate its influence on the result an agreement of 10 decimal digits was found between 
the numerical program and the analytical results (Eq. (|22|) and Eq.fl2"3"|)). 



3 Conclusions. 

An extensive use of the differential equation method allowed to find closed analytical 
expressions for the sunrise master integrals with one mass equal to zero, two arbitrary 
masses and an arbitrary external four momentum. 

Acknowledgments. 

One of us (HC) is grateful to M. Caffo and E. Remiddi for valuable discussions con- 
cerning the sunrise graph. 



References 

[1] F.V. Tkachov, Phys. Lett.B 100 (1981) 65; K.G. Chetyrkin and F.V. Tkachov, Nucl. 
Phys.B 192 (1981) 159. 

[2] F.A. Berends, M. Bdhm, M. Buza and R. Scharf, Z. Phys.C 63 (1994) 227. 

[3] M. Caffo, H. Czyz, S. Laporta and E. Remiddi, Nuovo Cim.A 111 (1998) 365, [Eep= 

| th/9805TT8; . 

[4] F.A. Berends, A.I.Davydychev, N.I. Ussyukina, Phys. Lett.B 426 (1998) 95, [hep 1 

| ph/97122Qc . 



9 



[5] 
[6] 
[7] 
[8] 
[9] 



M. Caffo, H. Czyz and E. Remiddi, Nucl. Phys.B 581 (2000) 274, |hep-ph/ 99 12501 
A.I. Davydychev and V.A. Smirnov, Nucl. Phys.B 554 (1999) 391, |hep-ph/ 9903328 
S. Groote and A.A. Pivovarov, Nucl. Phys. B 580 (2000) 459, |hep-ph/0003TT5 



M. Caffo, H. Czyz and E. Remiddi, Nucl. Phys.B 611 (2001) 503, |hep-ph/01030T4 
J. Fleisher, M.Yu. Kalmykov and A.V. Kotikov, Phys. Lett.B 462 (1999) 169, Kep- 



ph/990524S . 



[10] 

[11] 
[121 



M. Argeri, P. Mastrolia and E. Remiddi, |hep-ph/02021~23 . 

E.L. Ince, Ordinary differential equations, Dover Publications, New York, 1956. 

M. Caffo, H. Czyz and E. Remiddi, |hep-ph/0203256 . 



10 



